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ABSTRACT 

We present new Spitzer 3.6/ini observations of a sample of disk galaxies spanning over 10 magnitudes 
in luminosity and ranging in gas fraction from ~ 10% to over 90%. We use these data to test 
population synthesis prescriptions for computing stellar mass. Many commonly employed models 
fail to provide self-consistent stellar masses in the sense that the stellar mass estimated from the 
optical luminosity typically exceeds that estimated from the near-infrared (NIR) luminosity. This 
problem is present in models both with and without TP-AGB stars, but is more severe in the former. 
Self-consistency can be achieved if NIR mass-to-light ratios are approximately constant with a mean 
value near 0.5 Mq/Lq at 3.6/im. We use the Baryonic Tully-Fisher relation calibrated by gas rich 
galaxies to provide an independent estimate of the color-mass-to-light ratio relation. This approach 
also suggests that the typical 3.6^m mass-to-light ratio is 0.5 Mq/Lq (0.65 Mq/Lq in the if^-band) 
for rotationally supported galaxies. These values are consistent with a Kroupa IMF. 
Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: photometry 
— galaxies: kinematics and dynamics — galaxies: stellar content 



1. INTRODUCTION 

One of the most fundamental properties of a galaxy 
is its luminosity and the mass of the stars that pro- 
duce it. Our understanding of stellar evolution is suf- 
ficiently advanced that it should be possible to compute 
the luminosity produced by a stellar population ab ini- 
tio. Indeed, there exist in the literature various prescrip- 
tions for estimating stellar mass fr om observed colors or 
spectral energy distributions (e.g..lBell fc de Jons 



Bell et al.ll2003l:lPortinari et al.ll2004l: iZibetti etal 



2001 



2009; 



Into fc Portinarill2013H . 

A general expectation of population synthesis models 
is that the relation between mass and light is more nearly 
constant in the near-infrared (NIR) than in the optical 
part of the spectrum. This follows from basic consider- 
ations: recent star formation populates the upper main 
sequence with luminous, blue stars. These stars pro- 
duce copious amounts of optical light from little mass 
and lead short lives, leading to substantial perturbations 
on the average mass-to- light ratio T* of a galaxy. The 
degree to which this occurs depends on the IMF and the 
intensity of a star forming event relative to the total stel- 
lar mass already present. These effects combine to make 
the prediction of any particular galaxy's optical mass-to- 
light ratio uncertain at the factor of a few level. Young 
stars contribute rather less in the NIR part of the spec- 
trum, so one expects a closer relation between light and 
mass at these wavelengths. Empirically, the scatter in 
the Tully-Fisher relation de clines as one goe s from blue 
to red to NIR wavelengths (IVerheijenl[2001f ). consistent 
with the expected decrease in scatter in T*. 



In this paper we present new 3.6/im Sptizer Space Tele- 
scope photometry of a sample of galaxies spanning a large 
range (ten magnitudes) in luminosity. This wavelength is 
expected to be close to the sweet spot of nearly constant 
T*. We combine these data with optical colors and lu- 
minosities to check the predictions of several population 
synthesis models. 

We describe the data in |j2l In f}3]we compute stellar 
masses with various population synthesis models, and es- 
timate the mean NIR mass-to-light ratio that is required 
for each model to be self-consistent. In g] we utilize 
the Baryonic Tully-Fisher Relation (BTFR) to obtain 
an independent, kinematic estimate of stellar mass. We 
discuss our results in ^jH including some of the implica- 
tions for galaxy gas fractions in §5.11 We summarize our 
conclusions in ^ 

2. DATA 
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The SINGS ('Kmnic utt et al.l 1 2003': 'Dale et al.' "2001 
and THINGS (Walter et al.l 12008; de Blok et al.. .20081) 
surveys have demonstrated the utility of Sptizer 3.6/Ltm 
data for constraining the stellar components of star form- 
ing galaxies. Here we wish to sample disk galaxies over 
as large a range of physical properties as possi ble. To 
this en d, we combine the THINGS data of lde Blok et al.l 
(|2008f ) with new data from from two Spitzer programs. 
One cycle 5 project targeted galaxies with high quality 
rotation curves selected to increase the sampling of circu- 
lar velocities to both higher and lower mass than already 
present in THINGS. A cycle 7 snapshot program pro- 
vides additional photometry for low surface brightness 
galaxies. The combined sample spans ten magnitudes in 
[3.6] luminosity. All galaxies have at least some ancillary 
optical photometry, and a large subset also has kinematic 
information that can be used for the Tully-Fisher rela- 
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Table 1 

Galaxy Data 



Galaxy 


D 


Mv 




B 


- V 


V 


- / 


J 


- K 


Refs. 


D564-8 


6 


51 


— 12 


64 


— 15 


14 









93 






1,3 


KK98 251 


6 


08 






— 16 


39 














1 


DDO 154 


4 


04 


— 14 


55 


— 16 


41 





32 









34 


2,3 


D631-7 


5 


49 


— 14 


50 


— 16 


73 





41 





55 






1.4 


D640-13 


13 


44 


— 14 


36 


— 16 


77 









55 






li4 


D568-2 


21 


3 


— 14 


6 


—16 


81 





45 





70 






1,3 


D572-5 


14 


56 


— 14 


56 


—16 


82 





44 





52 






li4 


F415-3 


10 


4 


— 15 


2 


—16 


99 





62 










1,5 


D500-3 


18 


52 


— 15 


74 


— 17 


18 









31 






1>4 


DDO 168 


4 


25 


— 15 


70 


-17 


45 





32 










1,3 


D575-1 


9 


79 


— 14 


22 


— 17 


51 









70 






1.4 


D512-2 


14 


12 


— 15 


22 


— 17 


83 









80 






1,4 


F611-1 


25 


5 


— 15 


4 


— 17 


92 





57 










1,5 


D500-2 


17 


89 


— 16 


38 


— 18 


25 





52 





42 






1,4 


DDO 170 


16 


3 






— 18 


50 














1 


F565-V2 


55 


1 


— 16 


2 


— 18 


76 





44 










1,5 


NGC 2366 


3 


27 


-17 


71 


— 18 


90 





54 









84 


2 


D723-5 


27 


7 


— 16 


9 


-19 


21 





55 





75 






1,3 


F563-V1 


57 


6 


— 17 


2 


— 19 


79 





23 





83 






1,5 


IG 2574 


3 


91 


— 18 


53 


—20 


15 





42 









58 


2 


F563-1 


52 


2 


— 17 


8 


—20 


40 





40 





86 






1,5 


F574-2 


92 


3 


— 18 


3 


—20 


50 





58 










1,5 


NGC 2976 


3 


58 


— 18 


33 


—20 


52 





55 









79 


2 


F568-V1 


84 


8 


— 18 


6 


—20 


82 





47 










1,5 


F561-1 


69 


8 


— 18 


3 


—20 


88 





69 










1,5 


F577-V1 


113 




— 18 


7 


—20 


95 





50 










1,5 


NGC 1003 


10 


24 


—20 


38 


—21 


12 





42 









73 


2 


UGC 5005 


57 


1 


-18 


7 


—21 


18 





35 










1,5 


F574-1 


100 




— 19 


1 


—21 


37 





51 










1,5 


F568-1 


95 


5 


— 18 


9 


—21 


38 





52 










1,5 


NGC 7793 


3 


61 


— 19 


42 


—21 


46 





63 









68 


2 


UGC 128 


58 


5 


— 19 


3 


—21 


88 





63 





68 






1,5 


NGC 2403 


3 


16 


— 19 


82 


—21 


97 





39 









75 


2 


NGC 925 


9 


43 


—20 


54 


—22 


30 





50 









83 


2 


NGC 2903 


8 


9 


—21 


48 


-22 


74 





55 









91 


2 


NGC 3198 


13 


8 


—21 


18 


—23 


00 





43 









92 


2 


NGC 3621 


6 


56 


-20 


57 


-23 


04 





52 









83 


2 


NGC 3521 


8 





-21 


62 


-24 


19 





68 









93 


2 


NGC 3031 


3 


65 


-21 


55 


-24 


28 





82 









88 


2 


NGC 5055 


8 


99 


-21 


76 


-24 


60 





64 









95 


2 


NGC 2998 


68 


3 


-22 


36 


-24 


77 





45 









99 


2 


NGC 2841 


14 


1 


-21 


95 


-24 


88 





74 









93 


2 


UGC 2953 


15 


4 






-25 


08 













90 


1 


NGC 6674 


51 


9 


-22 


17 


-25 


17 





57 









86 


2 


NGC 7331 


14 


9 


-22 


30 


-25 


30 





63 






1 


00 


2 


NGC 801 


75 


3 


-22 


30 


-25 


33 





61 






1 


05 


2 


NGC 5533 


59 


4 


-22 


16 


-25 


47 





77 









94 


2 


UGC11455 


75 


4 






-25 


61 










1 


16 


1 


UGCA 145 


64 


5 






-25 


61 










1 


10 


1 


UGC 2885 


75 


9 


-23 


30 


-25 


94 





47 









88 


1 



References. — Spitzer [3.6] magn itudes: 1 . Thi s w ork 2 . de Blok.e^.a]J 
J2008). Ancilla ry data: 3. .S chombc rFet al.l 1)20111 ) 4. nVachtcrnach et al] 
J2009I) 5. lMcGaueh fc de BloM 81998) . 

Note. — Galaxy photometric data in order of increasing [3.6] luminosity. 
Adopted distances are in Mpc. Hq = 75 kms~^ Mpc~^ is assumed when no 
direct determination is available. Absolute magnitudes assume = 4.83 

and M^-®' = 3.24. 



tion. 

The [3.6] luminosities of galaxies in the THINGS 
sample have been adopted from the mass models of 
Ide Blok et al.1 (|2008[ ). Only the total luminosity is used 
here; no distinction is made between bulge and disk com- 
ponents. For the remainder of the assembled sample, the 
Spitzer data have been analyz ed with the ARCHANGEL 
surface photometry package (|Schomberti 2011). Ellipti- 
cal isophotes have been fit to the data and integrated 
magnitudes determined from asymptotic fits to curves 
of growth. Special care has been taken to exclude fore- 



ground stars and replace the masked region with an es- 
timate of the galaxy light based on surrounding pixels. 
IRAC is a sensitive instrument, and many foreground 
stars (presumably copious numbers of M dwarfs) are de- 
tected that are not obvious in survey material. Careful 
cleaning of these contaminants is essential to accurate 
photometry. 

The data are presented in Table [1] For each galaxy, 
we assign a distance and compute the corresponding ab- 
solute magnitude in the optical y-band and NIR Sptizer 
IRAC [3.6] band. Distances are taken from the direct 
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Table 2 

Population Synthesis Mass-to-Light Ratios 



Model 


IMF 


av 


bv 






■Y-V 
^0.6 




-^[3.6] 
-■-0.6 


TP 


Bell & de Jong (2001) 


Scaled Salpeter 


-0.734 


1.404 


-0.692 


0.652 


1.28 


0.50 


0.41 


no 


Bell et al. (2003) 


Scaled Salpeter 


-0.628 


1.305 


-0.206 


0.135 


1.43 


0.73 


0.62 


no 


Portinari et al. (2004) 


Kroupa (2002) 


-0.654 


1.290 


-0.736 


0.730 


1.32 


0.50 


0.41 


no 


Zibetti et al. (2009) 


Chabrier (2003) 


-1.075 


1.837 


-1.390 


1.176 


1.07 


0.21 


0.14 


yes 


Into & Portinari (2013) 


Kroupa (2002) 


-0.900 


1.627 


-1.020 


1.054 


1.19 


0.41 


0.33 


yes 



Note. — Stellar mass-to-light ratios in the V and X-bands from various population synthesis models are 
given in solar units by the formula log = ai + bi{B — V). For reference, the mass-to- light ratios predicted 
by each model for B ~V = 0.6 are given. The last column denotes whether or not the model includes tracks 
for TP-AGB stars. 



measurem ents tabulated in the Extragalactic Distance 
Database (jTuUv et alj 120091 ) when available. When no 
direct distance measurement is known, a Hubble flow 
distance assuming Hq = 75 kms^^Mpc^^ is adopted. 
Measured colors are given where known. The majority 
of the sample have observed B — V colors. For brighter 
;alaxies, J — K ca n readily be extracted from 2MASS 
Jarrett et al.ll200dl) . The lower surface brightness galax- 
ies are typically not detected by 2MASS. Fo r many, V — I 
has been observed (jSchombert et al.ll2011[ ). References 
to the sources of these data as well as gas masses and 
kinematic data are given in the final column of Table [1] 



3. STELLAR MASSES FROM POPULATION SYNTHESIS 
MODELS 

Our knowledge of stellar evolution is sufficiently ad- 
vanced to enable the ab initio calculation of the spectral 
energy distribution (SED) of complex stellar populations. 
Considerable effort has gone into the development of stel- 
lar population synthesis models that do just this (e.g., 
iBruzual fc CharlotI [200l lMarastonll2005D . Outstanding 
progress has been made, and it has become standard 
practice to quote stellar masses for galaxies based on fits 
to multi-color data with such models. 

The accuracy with which stellar masses can be esti- 
mated is debatable, but is probably no better than a fac- 
tor of two. Outstanding problems include uncertainty in 
the IMF, variations in the star formation history of galax- 
ies, and the contribution of stars in brig ht but short-lived 
phases of evolution (e.g., TP-AGB stars : [Maras ton 2005). 
Here we compare the predictions of various models with 
each other. We also check their internal self-consistency 
to gage the extent to which the same model predicts the 
same stellar mass for the same galaxy when the luminos- 
ity is measured in different bands. 

A simple approach to estimating the stellar mass of 
a galaxy is to assume a single, constant mass-to-light 
ratio such that M^, = T^L. This is a crude ap- 
proximation, as we expect the mass-to-light ratio of a 
population to vary with age and, to a lesser extent, 
with m etallicity. For examp l e, usi ng a mutli-metallicity 
model (jSchombert fc RakosI [20091) . we find that a 12 
Gyr old stellar population of solar metallicity has Tjf = 
2.8 M0/L0, while a stellar population with the same age 
but peak [Fe/H] = -1.5 has Tjf = 1.8 Mq/Lq. That 
same solar metallicity stellar population has of only 
0.4 M0/L0 at an age of 1 Gyr. 

Hence a more sophisticated approach is to use a pop- 
ulation synthesis model to relate the mean mass-to-light 



ratio in band i to a color through 

log — ai + 6iColor. 



(1) 



This reduces to the simple approach if the slope b is small 
for some band i, as might be expected in the NIR. For the 
models mentioned above, we find that the solar metallic- 
ity model changes in color as it ages from 1 to 12 Gyr by 
A{B -V)= 0.37 and A( J - K) = 0.03. We thus expect 
more variation of T, with optical colors than with those 
in the NIR, at least for mature stellar populations. 

There should be some scatter about the mean T*-color 
relation. The scatte r in T<, is also expec ted to be min- 
imized in the NIR (jBeh fc de Jonell2001[ ). The scatter 
in the T^-color relation ultimately limits the accuracy 
achievable by this approach, even if all else if perfect. 
One might hop e to do better by using multi-color in- 
formation (e.g.. iZibetti et al1l2009| ). or fitting the entire 
SED. It is not obvious to what extent additional color 
information can reduce the scatter in T* given the de- 
generacy between age and metallicity. In any event, the 
accuracy of stellar masses inferred from SED fitting is 
limited by the quality of the population synthesis model 
that is employed. This approach will obviously be sub- 
ject to systematic errors if a model differs systematically 
from reality with wavelength. 

The coefficients and bi are given in Table [5] for sev- 
eral representative models for the i = V and K bands 
and B — V color. These particular choices are made be- 
cause the most data are available in these bands. The 
models typically stop at K while we now have a good 
deal of Spitzer [3.6] photometry. We relate the popula- 
tion synthesis predicted 7^-band mass-to-light ratio 
to T^^' using the model resuffs of lOh et al.l ()2008f ): 



T 



[3.6] 



0.92Tf 



0.05. 



(2) 



When it is necessary to make this conversion for the data, 
we assume Kg — [3.6] = 0.31. Thi s is the mean colo r that 
we obtain for the 74 galaxies of iDale et al.l (|2005f). We 
assume this value also applies to the many LSB galaxies 
in our data for which no high quality if-band data are 
available. Any mean shift in the Kg — [3.6] color between 
HSB and LSB populations is probably modest, and may 
well be masked by the scatt er, which is only 0.11 mag. 
in this color in the sample of IDale et al.l (|2005l ). 

Fig. [1] shows the V and [3.6] band mass-to-light ratios 
for the models given in Table [21 Rather than simply show 
a line for each model, we plot the actual galaxy data to 
emphasize the beads-on-a-string nature of this approach 
to estimating stellar mass: a single color may provide 
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Figure 1. The relation between B — V color and t lie ste llar mass-to-light ratio in the V-band (left) an d the Spitzer [3.6 ] band (right) from 
the population synthesis models of'Bell & de Jonp) (I2001J) with their scaled Salpetor IMF (open circles), IBell et al.' f2003) (filled circles; also 
scaled Salpctcr), Portinari ct al. (2004) with a Kroupa (2002) IMF (squares), Zibctti ct al. (2009) with a Chabricr (2003) IMF (triangles), 
and Int o &: Portinari (2013) with a , Kroupa (2002) IMF (stars). The formula of iOh et al.. (.2008. . equation [2} is used to convert between Ks 
and [3.6]. Note the large disparity between models in the NIR. 



a reasonable estimate of the mean mass-to-light ratio, 
but it cannot reproduce the scatter th at one expects. 
This is estimated to be from 0.1 dex (|Bell fc de Jongj 
I2001D to 0.15 dex (jPortinari et al.ll2004[ ) in the isT-band, 
and larger in the optical bands. Even this is an over- 
simplification, as the scatter may be a function of color, 
with larger scatter likely in very blue, actively star form- 
ing systems. Additionally, the mean color-mass-to-light 
ratio relation can bend (the slope becomes much steeper 
hv B -V < 0.55 in the mod els of iPortinari et aIll2004D 
or even bifurcate: the line of IBell et al.l ( 20b3f ) splits the 
difference between distinct branches of high and low 
at blue colors in their Fig. 20. 

Different models will differ if they adopt different 
evolutionary tracks or a different IMF. Both matter 
for the models consi dered here. As pointed out by 
IBell fc de Jong! (|200H ). changes in the IMF to include 
more or fewer low mass stars serve mostly to change the 
mass without much altering the luminosity or color, so 
to a decent approximation can be treated as shifts in 
T^: For specifi city, we adopt the scaled Salpeter IMF 
of IBell fc de Jo ng ( 200fl) for their model and also for 
that o f Bell ct jl| (|2003f ). We adopt the Kroupa IMF 
(jKroup a 2002 ) for the mod el of Portinari et aL ( 2004 ) 
and ISito fc Portinaril (120131). while IZibetti et all (|2009[ ) 



uses the Chabrier IMF (|Chabrieil I2003D . These models 
all give a similar run of T* with color in the optical (left 
panel of Fig. [Ij, with small offsets for the different IMFs 
as well as other detailed differences. The agreement is 
less good in the NIR (right panel of Fig. [H . There is a 
huge disparity between the model of IBell et al.l ()2003[ ). 
which has a relatively flat slope bk with a high normal- 
ization tJ,^'^' > 0.55 M0/L0, and that of IZibetti et al.l 
(jlOOS:), which has a steep dependence on color (even in 
the NIR) and a low normalization T ^f '^^ < 0.3 Mo/L©. 
The model of llnto fc Portina ri (2013") is intermediate be- 
tween that of IZib etti et al. (,2009') and t he older models. 
The models of Beh fc de Jonel (|2001( ) and IPortinari et al.l 
([2OM) largely overlap in between these two extremes, 
but it is unclear whether we should take comfort in 
this apparent a g reeme nt as the two have different IMFs. 
IBell fc de Jon^ (|200lD give a 0.15 dex offset between 



their scaled Salpeter IMF and the Kroupa IMF, in the 
sense that the latter is lighter. However , the Kroupa cas e 
appears not to include brown dwarfs ()Bell et al.ll2003D . 
which would increase the mass by about 0.15 dex. Here 
we simply adopt the various published coefficients with- 
out trying to reconcile definitions of the IMF. 

The chief difference in evolutionary tracks between the 
various model s considered here is the inclusion of TP- 
AGB stars bv IZibetti erall (I2009D a nd Unto fc Portinaril 
using the models o f lMarastonI (I2005D . These stars 
are in the latest stages of evolution, being short-lived 
and rare, but quite bright. Their contribution to the 
integrated luminosity of stellar populations is most pro- 
nounced in the NIR, where they greatly enhance the pre- 
dicted luminosity of galaxies while doing little to alter the 
predictions of previous generations of models in the opti- 
cal portion of the spectrum. Thi s results in th e low m ass- 
to-light ra tios of the modeltQ of IZibetti et al.l (|200il l and 
[Into fc Portinari (2013) in the right hand panel of Fig.[TJ 

We use the models in Table [5] with the data in Ta- 
ble [T] to compute the stellar masses of sample galaxies. 
These are reported in Table [3] For each galaxy, we use 
the observed color to predict the mass-to-light ratio sep- 
arately in both the optical and NIR for each model. We 
then use the corresponding luminos ity to obtain astel- 
lar mass estimate. The m odels of B ell fc de Jond (|200lD 
and IPortinari et all ()2004D are barely distinguishable, so 
we tabulate only the latter. This results in eight distinct 
mass estimates for each galaxy: two for each of the four 
models. For a few galaxies where only V — I is avail- 
able and not B — V, we utilize the available color with 
the same parame ters a s the other m odels, when possible 
()Bell et al.l (|2003D and IZibetti et al.l f200a) do not tabu- 
late coefficients for V — I). 

We can now compare stellar mass estimate from differ- 
ent population synthesis models. We can also check each 
model for internal self-consistency. If all is well, the stel- 
lar mass estimated for the same galaxy will be the same 



^ IZibetti et al.l {20091) advocate using multiband colors to esti- 
mate mass-to-light ratios, and their single-color coefficients that are 
reproduced in Table [2] are only approximations made to facilitate 
the sort of comparison made here. 
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Table 3 

Stellar Masses from Population Synthesis Models 



BOS 



P04 



Z09 



IP13 



„ 1 




^[3-6] 


M^ 




M^ 


p^[3.6] 




mS?-*^' 


D564-8 






7.14 


6.99 






7.08 


6.86 


DDO 154 


7.54 


7.63 


7.51 


7.24 


7.27 


6.45 


7.37 


7.01 


D631-7 


7.64 


7.77 


7.61 


7.45 


7.41 


6.80 


7.50 


7.26 


D640-13 






7.20 


7.21 






6.95 


6.69 


D568-2 


7.73 


7.81 


7.70 


7.51 


7.52 


6.90 


7.60 


7.35 


D572-5 


7.70 


7.81 


7.67 


7.51 


7.49 


6.89 


7.57 


7.34 


F415-3 


8.19 


7.90 


8.16 


7.72 


8.08 


7.27 


8.12 


7.63 


D500-3 






7.35 


7.03 










DDO 168 


8.00 


8.04 


7.97 


7.66 


7.72 


6.87 


7.83 


7.43 


D575-1 






7.39 


7.68 






7.21 


7.37 


D512-2 






7.96 


7.92 






7.83 


7.70 


F611-1 


8.21 


8.27 


8.17 


8.06 


8.06 


7.56 


8.12 


7.94 


D500-2 


8.53 


8.39 


8.50 


8.15 


8.36 


7.61 


8.43 


8.01 


F565-V2 


8.36 


8.58 


8.33 


8.28 


8.15 


7.67 


8.23 


8.11 


NGC 2366 


9.09 


8.65 


9.06 


8.42 


8.93 


7.91 


9.00 


8.30 


D723-5 


8.78 


8.78 


8.75 


8.56 


8.63 


8.05 


8.69 


8.43 


F563-V1 


8.48 


8.96 


8.45 


8.51 


8.16 


7.51 


8.29 


8.23 


IC 2574 


9.26 


9.14 


9.23 


8.82 


9.04 


8.18 


9.13 


8.64 


F563-1 


8.95 


9.24 


8.91 


8.91 


8.71 


8.24 


8.80 


8.72 


F574-2 


9.38 


9.30 


9.35 


9.10 


9.24 


8.61 


9.30 


8.99 


NGC 2976 


9.35 


9.30 


9.32 


9.08 


9.20 


8.57 


9.26 


8.96 


F568-V1 


9.36 


9.41 


9.32 


9.14 


9.16 


8.55 


9.22 


8.98 


F561-1 


9.52 


9.47 


9.49 


9.34 


9.44 


8.94 


9.47 


9.27 


F577-V1 


9.44 


9.47 


9.40 


9.21 


9.26 


8.66 


9.33 


9.07 


NGC 1003 


10.00 


9.52 


9.97 


9.21 


9.78 


8.57 


9.87 


9.03 


UGC 5005 


9.24 


9.54 


9.21 


9.18 


8.98 


8.44 


9.08 


8.96 


F574-1 


9.61 


9.67 


9.58 


9.39 


9.43 


8.84 


9.50 


9.25 


F568-1 


9.54 


9.65 


9.51 


9.40 


9.37 


8.86 


9.44 


9.27 


NGC 7793 


9.89 


9.69 


9.86 


9.52 


9.78 


9.08 


9.83 


9.43 


UGC 128 


9.85 


9.86 


9.81 


9.69 


9.73 


9.24 


9.78 


9.60 


NGC 2403 


9.74 


9.86 


9.71 


9.53 


9.50 


8.85 


9.59 


9.33 


NGC 925 


10.17 


10.01 


10.14 


9.75 


9.99 


9.20 


10.06 


9.61 


NGC 2903 


10.61 


10.19 


10.58 


9.97 


10.46 


9.46 


10.52 


9.85 


NGC 3198 


10.34 


10.28 


10.30 


9.97 


10.12 


9.34 


10.20 


9.80 


NGC 3621 


10.21 


10.30 


10.18 


10.06 


10.04 


9.53 


10.11 


9.93 


NGC 3521 


10.84 


10.79 


10.80 


10.66 


10.75 


10.24 


10.79 


10.58 


NGC 3031 


10.99 


10.85 


10.96 


10.80 


10.98 


10.48 


10.99 


10.78 


NGC 5055 


10.84 


10.95 


10.81 


10.79 


10.74 


10.35 


10.78 


10.70 


NGC 2998 


10.83 


10.99 


10.80 


10.70 


10.63 


10.09 


10.71 


10.53 


NGC 2841 


11.05 


11.07 


11.01 


10.98 


11.00 


10.61 


11.02 


10.93 


NGC 6674 


10.92 


11.17 


10.88 


10.96 


10.77 


10.47 


10.83 


10.85 


NGC 7331 


11.05 


11.23 


11.01 


11.06 


10.93 


10.61 


10.98 


10.97 


NGC 801 


11.02 


11.24 


10.99 


11.05 


10.90 


10.59 


10.95 


10.96 


NGC 5533 


11.17 


11.32 


11.14 


11.24 


11.14 


10.89 


11.15 


11.20 


UGC 2885 


11.24 


11.46 


11.20 


11.18 


11.04 


10.60 


11.11 


11.03 


Note. — 


Masses are base ten logarithms in 


Mq. For each population model, 



the mass estimate from the V-band luminosity is given first, the n that from the 
[3.6] luminosity. The model s used are tho se of Bell ct al. ( 2003iV [Portinari et al.l 
1120041) . IZibetti et al.l | |2009| ). and Unto fc Portinari 12013) . 



irrespective of whether the luminosity is measured in the 
optical or NIR. Scatter in the color-mass-to-light rela- 
tion precludes this from ever being exactly true, but the 
sample is large enough that we can check whether it is 
true on average. This is done in Fig. [2l which shows the 
stellar mass estimated from the [3.6] luminosity plotted 
against that estimated from the y-band luminosity. 

The external consistency of the models is fairly good 
in the optical: examination of Table [3] shows that 
MY is usually similar across the board, with small off- 
sets that are roughly consistent with differences in the 
adopted IMFs and the intrinsic uncertainties of com- 
plex stellar populations. The same cannot be said in the 
NIR. Indeed, most of the models are not internally self- 
consistent: for a given stellar mass, they over-predict the 
infrared luminosity relative to the optical luminosity. In 
other words, the predicted infrared mass-to-light ratios 



are too small relative to those in the optical. 

The data cover many decades in luminosity. Conse- 
quently, even a small offset from the line of equality in the 
logarithmic Fig. [5] corresponds to a seriou s misestimation 
of the mass-to-light ratio. The models of Bell &: de Jon^ 
llOOlj) and iPort inari ct al. (2004) suffer this problem, 
but it becom e s mor e sev ere in the more recent m odels of 
iZibetti et all (I2009D a nd Unto fc Portinaril (|2013f ). In the 
case of iZibetti et all (|2009D . the offset in mass between 
optical and NIR estimates is typically greater than half 
a dex. 

The chief difference between the older and newer mod- 
els is the inclusion of TP-AGB stars in the more recent 
models. While it is obviously desirable to include all the 
relevant evolutionary phases, it appears that the con- 
tribution of the TP-AGB population to the integrated 
luminosity of galaxies is overstated by the models of 
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Figure 2. Stellar masses (Table [St e s timated from luminosities and colors using the formulae supplied by population synthesis models 
IIBell et al.l 12003 : IPortinari e'ral|[200l : IZibetti et alll2009l : Unto &: Portinaril [20Tj) . For each model, the mass estimated from the [3.6] 
luminosity of each galaxy is plotted against that estimated from the V-band luminosity. Filled circles use B — V colors to estimate mass- 
to-light ratios; open circles use V — I colors. If the models were perfect the data would follow t he line of unity, m odulo the expected 



scatter in the relation between the mass-to-light ratio and color. With the exception of the model of lBell et al.l ||2603I ) . the models tend to 
over-predict the NIR luminosity relative to the optical luminosity for a given stellar mass. 

between optical and NIR luminosity identified here will 
result in a systematic skew towards lower stellar masses 
as more NIR data are incorporated into SED fits. 
The model of Bel l et al.l (j2003f ) is not entirely ab initio, 



iMarastonl (I2005D. This conc lusion is consistent with the 
finding of lKriek et al.l ()2010D , who attempted to fit these 
models to the SEDs of post-starburst galaxies. Similarly, 
IZibetti et al.l (|2013f ) do not observe the spectral features 
expected from TP- AGB stars in the NIR sp ectra of post- 
starburst galaxies. IMelbourne et al.l (j2012D also reach a 
similar conclusion about the ratio of NIR to optical lu- 
minosity from the color-magnitude diagrams of nearby 
galaxies. 

Our understanding of the latest phases of stellar evo- 
lution does not yet appear to be sufficiently complete 
to confidently predict the NIR spectra of complex stel- 
lar populations. All ab initio models considered here 
suffer the same problem, which has gotten worse rather 
than better with good faith efforts to extend evolution- 
ary tracks to late times. Our results provide an empirical 
constraint to help inform future modeling. In the mean 
time, considerable caution is wa rranted in assigning ste l- 
lar masses based on SED fits (IConrov fc GunnI 1201(11 ). 
The results will depend not only on the model adopted, 
but also one what range of wavelength is fit. The offset 



as it makes extensive use of empirical constraints. This 
grounding in data makes it the only model to come close 
to the ideal of obtaining equal masses independently from 
optical and NIR luminosities. This might be a bit of a 
fluke given the bifurcation noted above, but nevertheless 
the ratio of optical to NIR fiux is approximately that 
observed in real galaxies. 

Referring back to Fig. [U we can immediately take two 
lessons: (i) self-consistency is only obtained for relatively 
high NIR mass-to-light ratios, and (ii) the slope of the 
color-mass-to-light ratio is not steep in the NIR. Con- 
sistency over the full range of the data only follows if 

xi^'^' is approximately constant. It has, of course, been 
suspected for a long time that T, would be roughly con- 
stant in the NIR. The data presented here support that 
expectation. 

We can perform a simple exercise to pursue the appar- 
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Figure 3. Stellar masses as in Fig. [2]but assuming a constant mass-to-light ratio at [3.6] rather than the population synthesis value. The 

values use d arc those that on average agree with the V -band stell ar mass: xj?'^^ = 0.58; 0.53; 0.43; and 0.48 Mq/Lq for the models of 
I Bell et al.l ||2003); Portinari et al. (2004); Zibetti ct aQ H2009I '): and Unto fc Portinaril (2013), respectively. Comparison to Fig. [2] makes it 
clear that a constant mass-to-light ratio in the NIR provides a more self-consistent description of the data than do models that apply a 
color correction in that part of the spectrum. 



ent constancy of xjf' . In Fig. [3] we retain the stellar 
mass predicted from the optical luminosity [M^ from 
Table[3]), but for the NIR based mass we adopt a constant 
mass-to- light ratio for all galaxies. We then ask what 

value of xi^'^' best reconciles the data with the optically 
estimated mass. The best fit mean NIR mass-to-light 

ratio is T^^' = 0.52 ; 0.58; 0.53; 0.43: and 0.48 M^^/I^ 
for the models of (Bell fc de Jo ng (2001); Bel l et alj 
(|2003[) ; |Portinari et al.l (|2004D ; iZibetti et al.i (|2009l) ; and 

Unto fc Portinaril ()2013[ ). respectively. Note that xi^'^' 
must be relatively large to obtain self-consistency even 
for the model with the lowest mean mass-to-light ratio 
(jZibetti et al.l l2009). 

Fig. [3] now looks like it should: a one to one correspon- 
dence between stellar mass estimated from optical and 
NIR luminosity, with the inevitable scatter in X* from 

galaxy to galaxy. Simply assuming a constant xlf '^^ per- 
forms better than the models consider ed here (compare 
Figures [2] and [31). Only in the case of iBeh et al] (|2003[ ) 
are the two not statistically distin guishab l e. As noted 
above, the relative success of the .Bell et al.l (pOOl model 



stems in part from its relatively flat slope, — 0.135. 
All of the other models listed in Table[3]have hK > 0.65, 
and two have &if > 1. Such steep slopes prevent self- 
consistency. There is no strong evidence in the data for 
a non-zero slope in the color-mass-to-light relation at 
3.6/im. A small slope is of course allowed, but may well 
be lost in the scatter in X*. 

The absolute value of the required x[f '^^ does not vary 
much from model to model. Given that these models 
bracket the plausible range in IMFs, there is tentative 
reason to hope that the absolute value of x[f'^' esti- 
mated in this way is approximately correct. Apparently, 



X 



[3.6] 



0.5 Mq/Jjq to a good approximation, with a 



range of about 20% from the Chabrier IIVIF at the low 
end to scaled Salpeter at the high end. 

The absolute value of the mass-to- light ratio given here 
implicitly assumes that the modeling of the optical por- 
tion of the spectrum is essentially correct. Really the 

data only constrain the ratio of xlf '^^ : TY to be higher 
than most models indicate. One is free to adjust the nor- 
malization of both optical and NIR mass-to-light ratio so 
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Figure 4. Left panel: the stellar mass estimated from the Baryonic TuUy-Fisher Relation (Table |4]| as a function of V-band luminosity. 
Right panel: the BTFR stellar mass as a function of [3.6] luminosity. Downward pointing triangles represent 2a upper limits for galaxies 



with stellar masses that are formally negative. The lines are the lines of equality for = 1.3 and xjf '^^ 
good proxy for kinematic stellar mass, with somewhat less scatter in the NIR than in the optical. 



0.5 Mq/Lq. Luminosity is a 



long as the appropriate ratio is preserved. However, ob- 
taining a NIR mass-to-light ratio much outside the range 

0.4 < T^^l < 0.6 M0/L0 implie s an IMEB ou tside the 
bounds of recent determinations (|Kroupall2002f ). 

4. STELLAR MASSES FROM THE BARYONIC 
TULLY-FISHER RELATION 

An independent method to estimate stellar masses is 
provided by the Baryonic TuUy-Fisher Relation (BTFR) . 
The observed mass of baryons — stars plus gas — 
is observed to correlate well with ro tation velocity 



(IMcGaugh et al 



I2001t IMcGaueh 



2000t IVerheiienI [200l IBell fc de Jonel 

Recent w ork on gas domi- 



2005) 



nated galaxies (iBeeum et al.l 120081: iStark et al.l 120091: 
iTrachternach et al.l I2009D has made it possible to ob- 
tain an absolute calibration of the mass scale of the 



BTFR that is largely independent of st ellar mass esti- 
mates (jStark et al.ll2009l:IMcGaughll2011[) . This provides 
a powerful new constraint on stellar masses. 
The BTFR can be written as 



Mh = AVf 



(3) 



where Vf is the rotation velocity measured in the flat 
outer portion of a resolved rotation curve and the bary- 
onic mass Mb — + Mg. A recent calibration of 
the B TFR obtains v4 = 47 ± 6 Mg km"^ s'* (|McGaughl 
[20T1 ). The gas masses in this calibration were taken to 
be Mg — 1.33Mhi, so the same is done here for consis- 
tencyO 

^ Our experience with point source contamination in the Spitzer 
data confirms that there are indeed stars that only become no- 
ticeable in the NIR. However, these apparently are foreground M 
dwarfs, the profusion of which anecdotally suggests that the IMF 
is not lacking in low mass stars. 

^ A few of the sample gala xies have mo lecular gas masses esti- 
mated from CO observations I ILerov et al.ll200^ . We have experi- 
mented with including molecular gas from this and from estimates 
based on the observed star formation rate. Unfortunately, such 
data do not exist for all sample galaxies. For those that do, the re- 
sult is essentially unchanged if we allow for an appropriately higher 
calibration of the BTFR {A = 50±6 M© km""' s'': McGaugh 2005). 
The stellar masses of individual galaxies change slightly, but the 
mean and scatter in the mass-to-light ratios are indistinguishable. 



In order to estimate stellar masses from the BTFR, 
we use the measured rotation velocity and equation [3] to 
first estimate the baryonic mass0 The stellar mass then 
follows by subtracting the gas mass: 



M, 



Mb - Mg = AVf - \MMhi 



(4) 



The data necessary to make this computation are 
the rotation velocity and the mass of atomic hy- 
drogen. T h ese a re taken f rom vari ous sources: 
Ide Blok et al.l (Il996l): I McGaugh^ (2005); IBegum et al.l 

Swaters et al.' (2009) 



(I2008D: Ide Blok ctal. (20081): 



IStark et al.. (2009); Tracht ernach et al., ((2009) and com- 
piled in Table ID The most restrictive requirement on 
the data is the need for a resolved rotation curve; this 
does not exist for all galaxies in Table [U but we have suc- 
ceeded in collecting a sample that covers a large dynamic 
range in Vf as well as in luminosity. 

The results of the computation of equation 2] are given 
in Table |4] as stellar mass and the mass-to-light ratios 

and xli^'^^ These are shown in Fig. 01 The stel- 
lar mass estimated in this way from the BTFR corre- 
lates as expected with luminosity, especially -/j[3.6], even 
though equation[4]involves only rotation velocity and gas 
mass, and makes no reference to the observed luminos- 
ity. Furthermore, the mean normalization is consistent 

with xii^'^' — 0.5 Mq/Lq. This provides a kinematic 
constraint on that is entirely consistent with what we 
inferred in the previous section. Stellar mass estimates 
from the BTFR are consistent with but independent of 
those based on requiring self-consistency in population 
synthesis. 

It is satisfying that the two independent methods for 
estimating stellar mass, population synthesis and the 
BTFR, provide convergent results. It is worth noting at 
this juncture a few exceptions that prove the rule. Look- 
ing at equation m it is remarkable that this works as well 
as it does. Starting from Vf which is typically of order 



* The empirical BTFR involves only the observed baryonic mass 
and the rotation velocity, and does not directly involve the mass 
of dark matter even if it is indire ctly implica t ed in setting the 
amplitude of the velocity See iMcGaughl l|20Tll . [20Tl l for further 
discussion. 
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Table 4 

Stellar Masses from the BTFR 



Galaxy 


Vf 




M* 


M, 


* 




-^[3.6] 

* 


C'^[3.6] 


D564-8 


29 


0.0210 


0.0126 


0.0228 


1.30 


2.34 


0.56 


1.00 


KK98 251 


40 


0.123 


-0.0027 


0.0392 






-0.04 


0.55 


DDO 154 


53 


0.420 


-0.0410 


0.142 


-0.73 


2.51 


-0.56 


1.95 


D631-7 


53 


0.195 


0.172 


0.145 


3.20 


2.69 


1.76 


1.48 


D500-3 


45 


0.0874 


0.112 


0.107 


0.66 


0.63 


0.76 


0.73 


DDO 168 


53 


0.421 


-0.0394 


0.0866 


-0.24 


0.53 


-0.21 


0.46 


D512-2 


38 


0.0926 


0.0003 


0.0659 


0.0032 


0.63 


0.0012 


0.24 


D500-2 


68 


1.15 


-0.147 


0.415 


-0.48 


1.36 


-0.37 


1.04 


DDO 170 


65 


1.14 


-0.302 


0.149 






-0.61 


0.30 


NGC 2366 


60 


0.798 


-0.201 


0.419 


-0.19 


0.40 


-0.28 


0.58 


IC 2574 


77 


1.88 


-0.249 


0.434 


-0.11 


0.20 


-0.11 


0.19 


F563-1 


111 


4.65 


2.48 


2.72 


2.20 


2.41 


0.86 


0.95 


NGC 2976 


86 


0.178 


2.39 


0.531 


1.30 


0.29 


0.75 


0.17 


F568-V1 


124 


5.22 


5.89 


3.85 


2.50 


1.63 


1.39 


0.91 


NGC 1003 


114 


6.78 


1.01 


1.12 


0.084 


0.093 


0.18 


0.20 


UGC 5005 


99 


7.04 


-2.52 


1.91 


-0.98 


0.74 


-0.43 


0.33 


F574-1 


100 


5.46 


-0.766 


1.97 


-0.21 


0.53 


-0.11 


0.28 


F568-1 


116 


5.34 


3.16 


3.12 


1.02 


1.00 


0.45 


0.44 


NGC 7793 


110 


1.01 


5.97 


1.42 


1.18 


0.28 


0.78 


0.19 


UGC 128 


131 


8.66 


5.177 


4.58 


1.15 


1.02 


0.46 


0.41 


NGC 2403 


136 


3.35 


12.6 


2.77 


1.74 


0.38 


1.03 


0.23 


NGC 925 


114 


5.31 


2.52 


2.00 


0.18 


0.14 


0.15 


0.12 


NGC 2903 


188 


5.78 


52.3 


10.3 


1.56 


0.31 


2.10 


0.42 


NGC 3198 


149 


13.5 


9.57 


3.93 


0.38 


0.16 


0.31 


0.13 


NGC 3621 


152 


9.28 


16.0 


3.86 


1.11 


0.27 


0.49 


0.12 


NGC 3521 


192 


5.96 


57.8 


24.7 


1.52 


0.65 


0.61 


0.26 


NGC 3031 


199 


9.27 


64.4 


22.5 


1.80 


0.63 


0.63 


0.22 


NGC 5055 


181 


9.59 


41.0 


14.7 


0.95 


0.34 


0.30 


0.11 


NGC 2998 


212 


36.8 


57.5 


12.9 


0.77 


0.17 


0.36 


0.08 


NGC 2841 


279 


11.4 


272. 


58.8 


5.28 


1.14 


1.52 


0.33 


UGC 2953 


260 


9.28 


205. 


59.6 






0.97 


0.28 


NGC 6674 


241 


50.8 


106. 


22.8 


1.68 


0.36 


0.46 


0.10 


NGC 7331 


246 


12.5 


159. 


32.0 


2.24 


0.45 


0.61 


0.12 


NGC 801 


219 


48.8 


58.9 


15.3 


0.83 


0.22 


0.22 


0.06 


NGC 5533 


240 


31.2 


124. 


23.7 


1.98 


0.38 


0.41 


0.08 


UGC11455 


272 


24.9 


232 


42.2 






0.67 


0.12 


UGCA145 


300 


52.1 


329. 


54.8 






0.94 


0.16 


UGC 2885 


298 


45.8 


324. 


68.6 


1.82 


0.39 


0.69 


0.15 



Note. — Stellar masses and mas s-to-light ratios inferred from the Baryonic 
TuUy-Fisher relation l lMcGaughll2012l 'l. The rotation velocity is in kms~^. Masses 
are in units of 10^ Mq . The gas mass assumes Mg = 1.33M}jj . Mass-to- light ratios 
arc in solar units Mq/Lq. 



~ 100 kms~^, we first raise this number to the fourth 
power. We then multiply by an empirically determined 
constant of the same order. This gives a large number 
appropriate to the baryonic mass of a galaxy. We then 
subtract another large number, the gas mass. Anyone 
with experience in attempting this sort of calculation 
with astronomical data will immediately recognize the 
dangers inherent in such an approach. Raising a number 
to a large power invites the propagation of errors; sub- 
tracting one large number from another practically begs 
them to blow up. Remarkably, this rarely happens. This 
is a credit to both the quality of the data and the virtue 
of the method. 

It does sometimes happen that the stellar mass esti- 
mate from equation|3]is negative. This is especially prone 
to occur in gas rich galaxies where the gas mass accounts 
for nearly all of the baryonic mass. There are nine oc- 
currences of this among the 38 galaxies listed in Table |31 
The formal significance of these negative masses is < 2a 
in all cases, so there are no examples where the method 
yields a definitively unphysical result. Indeed, it would 
be worrisome if one did not obtain the occasional nega- 



tive mas^, especially given the large number of gas rich 
galaxies in this sample (about half have fg > 0.5). 

At the opposite extreme, maximum disk places an up- 
per limit on the allowed stellar mass: the inferred mass 
should not cause the rotation curve of the stellar com- 
ponent to exceed the observed rotation curve. Since 
the BTFR only unitizes the outer part of the rotation 
curve, it is in principle quite possible for this to occur. 
In this sample, it occurs exactly once: in the case of 
NGC 2841. This is hardly surprising, as NGC 2841 
is a well known outl ier from the Tully-Fisher relation 
(jBottema et al.ll2002f) . If it is a true outlier, then this 
method does not apply. It remains conceivable that its 
distance is wrong or some other systematic error is at 
fault. For example, an undiagnosed warp in the outer 
regions could cause an error in Vf. Regardless of what 
is really going on in this case, the mass returned by the 
BTFR is unphysical in that it exceeds maximum disk. 

^ We have attempted this exercise using line-widths instead of 
Vf from rotation curves. Line- widths fail by giving many negative 
stellar masses and some unreasonably large T,. This happens be- 
cause of the larg er intrinsic scatter of the line-width based BTFR 
IIMcGaughl 12012 ): good, resolved rotation curve are essential for 
estimating T*. 
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Figure 5. Mass-to-light ratios inferred from the baryonic Tully-Fisher relation in the V-band (top) and at [3.6] (bottom) as a function of 
B — V (left) and J — Ks (right) color. Galaxies for which the uncertainty in the stellar mass exceeds half of the stellar mass are shown as 
smaller symbols. Downward pointing triangles represent 2(t upper limits for galaxies with formally negative T*. NGC 28 41 is shown as an 
open circle as the inferred stellar mass exceeds its maximum disk value. For comparison, the stellar population models of lBell et al.l I I2003I , 
solid lines) and lPortinari et al.l II2004I . dashed lines) are shown in the left panels. The median mass-to-light ratios are = 1.09 Mq/Lq 



and T^'^l 
ratio is T? 



: 0.45 M0/L0 (dotted lines in right panels). For the average Ks — [3.6] = 0.31 color, the corresponding i^s-band mass-to-light 
= 0.58 Mq/Lq. Note that the majority of galaxies are within 2a of the median value in the NIR. 



The only other galaxy that appears troublesome is 
NGC 2903. Its stellar mass is not unphysical, but it 
is odd in having its BTFR-inferred NIR mass-to-light 

ratio larger than the optical value: tI'^'^' > ■ Popula- 
tion synthesis models do not anticipate that this should 

happen, even after the upward correction to t1^'^' dis- 
cussed in ij3l nor does it: all other galaxies but twcQ 
have BTFR-inferred Tp'^' < T]f. Looking at Tabled 
the stellar masses for NGC 2903 from optical and NIR 
are badly offset according to all four models, as if the 
NIR luminosity had been underestimated. Indeed, from 
a purely empirical perspective, NGC 2903 is an outlier 
from the infrared luminosity Tully-Fisher relation even 



6 D500-3 and NGC 1003 have tI,^ '^' ^ T^, indistinguishable 
within the errors. More generally, population models do antici- 
pate low and rapidly variable mass-to-light ratios for populations 
dominated by very young stars, but this situation stabilizes to 
-jpp.e] ^ yY £qj. jjjQg^ relevant ages (> 1 Gyr). 



though it is not from the optical version of the relation . 
The source of the photometry is Ide Blok et all ()2008D . 
who note problems in constructing the mass model for 
this galaxy, so we regard its [3.6] luminosity as dubious. 

This is a small list of exceptional cases for a method 
that should be so easily derailed. Just as the stellar mass 
might come out negative, if could also come out an order 
of magnitude or more too heavy. That this does not often 
happen is encouraging. 

Fig. [5] shows the mass-to-light ratios inferred from the 
BTFR as a function of color. The data show a trend 
TY increasing with B — V color that is consistent 
with the expectations of population models. This has al- 
ready been seen to be the case in the B-band (jStark et al.l 
120091 ): the present data are entirely consistent with that 
earlier result. There is a fair amount of scatter in T*, 
but this is also expected. The inferred optical mass-to- 
light ratios behave as they should from the perspective 
of population synthesis. 

There is no clear indication of a correlation of either 
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Figure 6. Mass-to-light ratios inferred from the baryonic TuUy-Fisher relation as per Fig. (5] but restricted to the more precise data with 

""M* < ^^'^f The mean mass-to-light ratios are Xjf = 1.3 Mq/Lq and xjf '^' = 0.5 Mq/Lq (dotted line in right panel). For the average 

Ks — [3.6] color, the corresponding Ks-hand mass-to-light ratio is T^" = 0.65 Mq/Lq. This is slightly heavier than seen for all the data 
in Fig. [5] as the more precise data tend to be for the brighter, more evolved spirals. The line in the left panel represents a fit to the data 
with the slope fixed to the value anticipated by population synthesis models (Table [SJ. 



mass-to-light ratio with J — A'^ color. This is not sur- 
prising, since this color should be more s ensitive to the 
metallicity of stars than their age (jBell &: de Jong 2000). 
This means that NIR colors are a poor indicator of T,. 
One is better off assuming constant xlf '^^ than attempt- 
ing to correct T» with a color term involving NIR colors. 

The NIR mass-to-light ratio is also insensitive to op- 
tical color. The data allow for a small slope fep.e], 
but are also consistent with zero slope. This is con- 
sistent with the results of where it was seen that 
small slopes were required for population synthesis mod- 
els to be internally self-consistent. To a reasonable ap- 
proximation, xi^'^' is constant. The median mass-to- 
light ratio is xlf'^' = 0.45 M0/L0. For the average 
Ks — [3.6] = 0.31 color, the corresponding mass-to-light 
ratio is Xf = = 0.58 Mq/Lq. 

The scatter in Fig. [S] is somewhat misleading because 
of the large uncertainty in stellar mass estimated from 
the BTFR for gas rich galaxies. In Fig. [S] we restrict 
the sample to those galaxies with formal errors in mass 
o-M./M* < 50%. We also exclude NGC 2841, whose 
BTFR mass is unphysical, and NGC 2903, whose [3.6] 
luminosity we consider untrustworthy (see above). This 
leaves a "clean" sub-sample of 14 galaxies. 

The scatter in Fig. [5] is reduced from that in Fig. [5l 
as expected by construction. The 1^-band mass to-light- 
ratios remain consistent with the color trend expected 
from population models, but there remains some scat- 
ter. The [3.6] mass-to-light ratio is effectively constant. 
There is no hint of a trend with color, and many in- 
dividual galaxies are consistent with the average value 
within the errors. The average mass-to-light ratio is 
slightly higher for the clean sub-sample, with a mean 

Xi^-^l = 0.5 Mq/Lq (equivalent to Xf = = 0.65 Mq/Lq). 
This small shift from 0.45 Mq/Lq for the entire sample is 
too small to consider significant. Nevertheless, it goes in 
an intriguing direction, hinting that gas rich galaxies may 
tend to have lower mass-to-light ratios, even in the NIR. 
It is not difficult to imagine that they have younger ages 
and/or are more metal poor than bright, well-developed 
spirals. 



The clean sub-sample gives us some handle on the 
intrinsic scatter in mass-to-light ratios. The scatter 

in xi'^'^^ is 0.23 dex while the average uncertainty is 
0.19 dex. This implies an intrinsic scatter of 0.13 dex. 
While the sample is too small and the uncertainty in 
the uncertainties too great to draw firm conclusions from 
this, this implied value of the intrinsic scatter is nicely 
brac keted by that anticipa ted in the NI R by the mod- 
els of [Bell fc de Jond (|2001l . 0.1 dex) and iPortinari et al.l 
(1200J, 0.15 dex). 

Table [5] summarizes our calibration of the color — mass- 
to-light ratio relation determined from the BTFR. This 
is entirely consistent with the expectation of population 
synthesis models once those are made self-consistent as 
discussed in [J3] We do not have enough data to constrain 
the slope of the relation, which we fix to a representative 
value from the population models in the optical bands. 
The B-band calibration is given along with that in the 
t/-band for completeness; see also lStark et"all (2009). In 
the NIR there is no indication of a systematic variation 
of X* with color, so we fix the slope to zero. The best 
fit to the NIR X* is a constant, differing slightly by sam- 
ple selection. The clean sub-sample gives a more pre- 
cise constraint, but may build in a small systematic bias: 
the uncertainty on the stellar mass from the BTFR is 
necessarily smaller for star dominated rather than gas 
dominated galaxies. However, the difference is not large 
(Xi^'^' — 0.5 versus 0.45 Mq/Lq), and we do not con- 
sider it to be meaningful. 

The basic value of the BTFR const rained mass- t o-ligh t 
ratio is meaningful. As discussed by iStark et "all ()2009D , 
it provides a constraint on the IMF. Comparing ay de- 
termined from the BTFR in Table [5] with that expected 
from population synthesis models in Table [2l it appears 
that the data are most consistent with a Kroupa IMF 
(lKroupall200l) . 

5. DISCUSSION 

For all intents and purposes, the data appear to be 
consistent with a nearly constant NIR mass-to-light ra- 



tio X 



[3.6] 



0.5 Mq/Lq. This value is inferred in two 
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Table 5 

Self-Consistent Stellar Mass-to-Liglit Ratios 



Sample as bs ay by aK, bK, a[3.6] ^'[3.6] '^E.a '^ae "^ae "^0.6 

All galaxies -0.94 1.7 -0.68 1.3 -0.24 -0.35 1.20 1.26 0.58 0.45 

Clean sub-sample -0.90 1.7 -0.64 1.3 -0.19 -0.30 1.32 1.38 0.65 0.50 

Note. — Stellar mass-to-light ratios (in solar units) log T* = ai + bi{B — V) tied to the gas rich galaxy 
calibration of the BTFR (McGaugh 2012). The slope is fixed so that only the zero point varies with sample. 
The clean sub-sample consists of galaxies with relatively precise stellar masses estimated from the BTFR with 
formal uncertainties om, < Implied mass-to-light ratios a.t B — V = 0.6 are given for filters i = B, V, 

Ks, and [3.6]. 



independent ways: (i) by requiring self-consistency in 
population synthesis models, and (ii) from the stellar 
mass required to obey the BTFR hxed by gas dominated 
galaxies. The equivalent Kg-hand mass-to-light ratio is 
Tf = 0.65 Mq/Lq assuming K, - [3.6] = 0.31 with so- 

lar absolute ma gnitudes Mq" = 3.27 and m|^-^' = 3.24 
(|Oh et al.ll20"08h . 

The apparent near-constancy of the stellar mass-to- 
light ratio in the NIR that we infer here has long been an- 
ticipated. This observation has important ramifications. 
Simply by observing the NIR luminosity of a galaxy, we 
obtain a good handle on its stellar mass — no further 
input is required. 

It is tempting to conclude not only that the NIR lu- 
minosity of a galaxy is a good indicator of its stellar 
mass, but also that we have constrained the conversion 

of light to mass to within - 20% of T^^' = 0.5 Mq/Lq. 
There is certainly real scatter within that range, and 
there are hints in the data that the range 0.4 - 0.5 may 
be more likely among gas rich disks than the range 0.5 - 
0.6. While it is easy to imagine that some galaxies with 
genuinely young stellar populations have even lower T*, 
this does not appear to be the general rule. 

Our chief hesitation in declaring that the problem of 
stellar mass is solved (besides historical precedence that 
such ass ertions are never warr anted) is that the DiskMass 
survey (|Bershadv et al.llMO,) finds masses that a ppear 
to be systematically lower ((Bershadv et al.l l20lH ) by a 
factor of two or even a bit more. There is no overlap in 
the samples, so this estimate is inferred from their esti- 
mate of the maximality of disks, which is small, whereas 
our T, are more in line with, or even a bit heavier than 
those of Ide Blok et al.l (|2008t see then Figure 61). We 
do not know the reason for this inconsistency. 

We can note some consequences for our study of a much 
lower mass-to-light ratio. From the stellar population 
perspective, it would be possible, at least in principle, 
to reduce T, of both optical and NIR bands while pre- 
serving the ratio found in fJ31 In practice this is difficult. 
The optical port ion of the spectru m is reasonably well 
understood fe.g. JKriek et al.ll2010D . so such a reduction 
would correspond to a large shift in the IMF. That is not 
inconceivable, but would have real consequences for star 
counts: we cannot reduce the IMF arbitrarily and still 
account for the stars we see directly. 

It is also difficult to reconcile lower with the BTFR. 
The observed slo pe of the i^-band TuUy-Fisher relation 
(|Verheiienll2001[ ) is 4, c onsistent with that of gas dom- 
inated galaxies (jMcGau gh 2012). The two join contin- 
uously if and only if the stellar mass-to-light ratio is as 



we estimate. If a lower mass-to-light ratio is adopted, 
the BTF R develops a brea k in its slope and the scatter 
mcreases (IMcGauehl 120051) . There is no hint of a break 
in the empirical relation, so it would be odd to insert 
one. It would also be odd if the scatter in the mass- 
based BTFR were larger than that in the luminosity- 
based TuUy-Fisher relation, since variations in T.^ should 
only add scatter when going from the intrinsic BTFR to 
the luminous Tully-Fisher relation. This would be es- 
pecially strange since the expected increase in scatter 
due to as one goes from the NIR t o progressively 
bluer bands is observed (jVerheiienll2001h . We therefore 
see no way to grac efully reconcile our results with those 
of iBershadv et al.l ([2011), so there appears to remain of 
systematic uncertainty in stellar mass of a factor of ^ 2. 

5.1. Gas Fractions of Disk Galaxies 

Gas fractions provide an i mportant clue to the evolu- 
tionary state of galaxie s fe.g. JMcGaugh fc de Blok 19971: 
iSchombert et al.l l2001[ ). As an example of the utility 

of a constant mass-to-light ratio, we adopt xi^'^' = 
0.5 M0/L0 in order to plot the gas fraction of galax- 
ies fg — Mg/{M^, + Mg) ^.s 3, fuuctlou of stellar mass 
in Fig. [T] This is particularly useful for gas rich galax- 
ies that have small but highly uncertain stellar masses. 
Fig. [7] shows a reasonable distribution of galaxies in the 
fg-M^ plane. Simply taking a constant NIR T, appears 
to give plausible stellar mass estimates, even for objects 
with formally negative stellar masses according to the 
BTFR method. 

Fig. [7] has a well defined upper envelope. This is easily 
anticipated from the definition of fg, which depends on 
the stellar mass. It is also real: galaxies with large stellar 
masses and high gas fraction are intrinsically rare. Such 
objects would be bright in both optical and 21 cm lumi- 
nosity, and readily detectable in large surveys. Objects 
with Af* 10^^ Mq have converted most of their origi- 
nal gas into stars by the present epoch. This is consistent 
with the downsizing of galaxies observed over a range 
of redshifts, in which massive galaxies typically appear 
to be the oldest and most evolved at any given epoch. 
Note, however, that this conclusion would not persist if 
were much smaller, as even the most massive spirals 
would then be inferred to be rather gas rich. 

Galaxies of low stellar mass run the gamut from very 
gas rich to gas poor. The smallest galaxies discussed 
here overlap with the classical dwarf Spheroidals of the 
Local Group in terms of stellar mass. Fig. [7] shows For- 
nax as an example of such objects. Though the lowest 
mass galaxies discussed here may have fewer stars than 
Fornax, they are still larger in terms of baryonic mass 
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Figure 7. The gas fraction as a function of stellar mass assuming 
~^13.6] _ Q J, jyjg^Lg The galaxy sample spans a large range in 
mass and samples a range of gas fractions at each mass. The ab- 
sence of galaxies with large stellar masses and large gas fractions 
is presumably real. Downsizing is evident in that the more mas- 
sive galaxies have typically consumed a larger fraction of their gas 
supply. One should however be cautious about the absence of low 
stellar mass galaxies with low gas fractions: this is presumably a 
selection effect, as HI must be detected to measure a gas fraction. 
Low mass but gas poor galaxies like the Fornax dwarf spheroidal 
lie in the lower left corner of this diagram. 

since their gas fractions are high. 

The paucity of galaxies in the lower left part of Fig. [7] is 
almost certainly a selection effect. In order to be included 
in our sample, a galaxy must have a positive HI detec- 
tion. Dwarfs like Fornax would not be included even if 
detected optically. Hence there is every reason to expect 
that galaxies fill in the space below the upper envelope. 
Indeed, our current sample does a good job of sampling 
the available parameter space, in spite of the inevitable 
cut-off at low masses and low gas fractions. 

The mere existence of gas rich disk galaxies with 
M, « 10^ Mq has interesting implications. This is not 
much larger than a globular cluster. It is widely ex- 
pected that the early re-ionization of the universe would 
suppress the accretion of gas into low mass halos (e.g., 
iRicotti fc GnedinI I2005D . This is frequently invoked to 
explain the properties of the dwarf spheroidals of the Lo- 
cal Group, in particular that they have no gas now, nor 
would have had for most of thei r history. If re-ionizatio n 
is the only part of this story (see lMcGaugh fc Wolj|2010[ ). 
it should hold everywhere: there should be no galaxies 
as low in mass as the dSph population with lots of gas. 
Yet there are examples in Figure[7]of objects with stellar 
masses c omparable to Fornax that are replete with gas 
(see also iTikhonov fc KlvpinI [200i ) . 

6. CONCLUSIONS 

We have presented new Spitzer data for disk galax- 
ies spanning ten magnitudes in [3.6] absolute magnitude. 
These galaxies also span the full range of gas fractions, 
from ^ 10% to over 90%. We have used these data to 
test the self-consistency of stellar population synthesis 
models from the optical to NIR bands. We have also 
applied the BTFR to obtain an independent, kinematic 
estimate of the calibration of stellar mass-to-light ratios. 

Our main conclusions can be summarized as follows: 

• Many commonly utilized stellar population models 



are not self-consistent in the sense that applica- 
tion of the same model to the same galaxy results 
in different stellar masses depending on whether 
an optical or NIR luminosity is used. This prob- 
lem occurs in the sense that the NIR luminosity is 
overestimated relative to the optical luminosity for 
a given mass of stars. This problem is present in 
all ab initio models considered, but is particularly 
severe in models that include TP-AGB stars. 

• Self-consistency between optical and NIR obser- 
vations can be achieved if NIR mass-to-light ra- 
tios are approximately constant and do not vary 
perceptibly with color. The typical value for self- 
consistency is 0.5 M0/LO at 3.6/^m (equivalent to 
0.65 M0/LO in the K^-band). 

• The BTFR calibrated by gas dominated galaxies 
provides an independent constraint on stellar mass. 
This approach also yields a typical 3.6/Ltm mass-to- 
light ratio of 0.5 Mq/Lq. 

• The intrinsic scatter about the typical mass-to- 
light ratio is modest for the most precise data, 
- 0.13 dex in the NIR. 

• The IMF is consistent with that of lKroupal ()2002D . 

Both independent methods yield the same typical mass- 
to-light ratio, giving some hope that this is the correct 
conversion between light and stellar mass. Within the 
context of our own study, it appears that the typical 
value is good to about 20%. We cannot exclude the pos- 
sibility of systematics at the factor of two level, though 
these appear unlikely. 

Workers wishing to estimate the stellar masses of 
galaxies would do well to adopt a constant NIR mass- 
to-light ratio as calibrated here. There appears to be 
little added value in fitting the complete SED. Indeed, 
such fits can only be as good as the population model 
the data are fit to. Given the systematic offsets in the 
models pointed out here, such SED fits are bound to suf- 
fer from systematic errors in stellar mass that will depend 
on the specific model employed and also on the range of 
wavelength fit. 
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the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under a contract with NASA. Support for this 
work was provided by NASA through an award issued 
by JPL/Caltech. Other aspects of this work were sup- 
ported in part by NASA ADAP grant NNX11AF89G and 
NSF grant AST 0908370. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California 
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